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1 Introduction 

Inflation is one of the leading paradigms describing the physical conditions that prevailed in 
the very early Universe [1—6]. It is a phase of accelerated expansion that solves the puzzles 
of the standard hot Big Bang model, and it provides a causal mechanism for generating 
scalar [7-11] and tensor [12] inhomogeneous perturbations on cosmological scales. These 
inhomogeneities result from the parametric amplification of the vacuum quantum fluctuations 
of the gravitational and matter fields during the accelerated expansion. 

The transition from these quantum fluctuations to classical but stochastic density per¬ 
turbations [13-16] gives rise to the stochastic inflation formalism [17-25]. 1 It consists of 
an effective theory for the long-wavelength parts of the quantum fields, which are “coarse 
grained” at a fixed physical scale (i.e. non-expanding), somewhat larger than the Hubble ra¬ 
dius during the whole inflationary period. 2 The non-conunutative parts of this coarse grained 
field tp are small, and at this scale, short-wavelength quantum fluctuations have negligible 
non-commutative parts too. In this framework, they act as a classical noise on the dynamics 

1 This formalism was, in fact, first used in Ref. [9] at the level of the Langevin equation, from which results 
lying beyond the one-loop approximation for the inflaton field were obtained. 

2 More precisely, the coarse grained part of the field consists of the modes k for which k < aaH. Here, a 
is a cutoff parameter satisfying [25] e -1 ^ 361 -* -C er -C 1, where t\ is the first slow-roll parameter. Under this 
condition, the physical results are independent of er. 
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of the super-Hubble scales, and p can thus be described by a stochastic classical theory, 
following the Langevin equation 


dp 

d N 


V' 
3 H 2 


+ i 


(i.i) 


This equation is valid at leading order in slow roll. Time is labeled by the number of e- 
folds N = In a, where a is the scale factor. The Hubble parameter H = da/{adt) is related 
to the potential V via the slow-roll Friedmann equation H 2 ~ V/(3M 2 l ), where M P1 is 
the reduced Planck mass. The dynamics of p is then driven by two terms. The first one, 
proportional to V’ (where a prime denotes a derivative with respect to the inflaton field), 
is the classical drift. The second one involves a white Gaussian noise, £, and renders the 
dynamics stochastic. It is such that (£ (IV)) = 0 and (£ (N) £ (N 1 )) = 6 (N — N'). 

The stochastic formalism thus accounts for the quantum modification of the super- 
Hubble scales dynamics. It allows us to calculate quantum corrections on background quan¬ 
tities beyond the one-loop approximation for the inflaton scalar field <j) (in fact, beyond 
any finite number of inflaton loops) and to calculate such quantities as e.g. the probability 
distribution and any moments of the number of inflationary e-folds in a given point. In 
turn, cosmological perturbations are affected too, and a natural question to address within 
the stochastic framework is therefore how quantum effects modify inflationary observable 
predictions. This is the main motivation of the present work. 

Stochastic inflation is a powerful tool for calculating correlation functions of quantum 
fields during inflation. In Refs. [25-29], it is shown that standard results of quantum field 
theory (QFT) are recovered by the stochastic formalism for test scalar fields on fixed in¬ 
flationary backgrounds, for any finite number of scalar loops and potentially beyond. This 
result is even extended to scalar electrodynamics during inflation in Refs. [30, 31] and to 
derivative interactions and constrained fields in Ref. [32]. In Ref. [26], fluctuations of a non¬ 
test inflaton field have been studied, too. In this last case, the calculation is performed at 
linear order in the noise, that is, assuming that the distribution of the coarse grained field 
remains peaked around its classical value <^ c , where </> c is the solution of Eq. (1.1) without 
the noise term. However, it may happen that the quantum kicks dominate over the classical 
drift and provide the main contribution to the inflationary dynamics in some flat parts of 
the potential. It is therefore legitimate to wonder what observable imprints could be left in 
such cases. In order to deal with observable quantities, the goal of this paper is therefore 
to calculate the correlation functions of inflationary perturbations in full generality, taking 
backreaction of created inflaton fluctuations on its background value into account, starting 
from Eq. (1.1) and without relying on a perturbative expansion in the noise. 3 

This work is organized as follows. In section 2, we first discuss the issue of the time 
variable choice in the Langevin equation (we further elaborate on this aspect in appendix A). 
This allows us to set a few notations, and to already argue why some of the effects later 
obtained (but not all) are Planck suppressed. In section 3, we turn to the calculation of the 
correlation functions of primordial cosmological perturbations, without assuming them to be 
small. We first review different methods that have been used in the literature, and motivate 

3 In this connection, the approach of Ref. [33] is close to ours. However, we use a different form of the 
Fokker-Planck equation, a different initial condition for the inflaton probability distribution, and a different 
form of the 5N formalism (which, in fact, may be called N formalism) that does not use an expansion in 5N 
and in the metric perturbation £ (in fact, these two quantities are not small in the so called regime of eternal 
inflation). 
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our choice of combining the stochastic and 5N formalisms. We then settle our computational 
strategy and proceed with the calculation itself. Results are presented in section 4; see in 
particular Eqs. (4.1) and (4.6). We show that the standard formulas are recovered in a 
“classical” limit that we carefully define, and discuss the regimes where they are not valid. 
Finally, in section 5, we summarize our main results and conclude. 


2 Time Variable Issue 


Because of the Friedmann equation, the Hubble parameter appearing in Eq. (1.1) is sourced 
by the inflaton field itself, through the slow-roll function H(ip). At leading order in the 
noise, one simply has H(ip) = H((j) c ), which is a classical quantity. Beyond the leading 
order, however, H is dependent on the full coarse grained field and is therefore a stochastic 
quantity. This has two consequences. The first one is that starting from a classical time 
label, any other time variable defined through a or H is a stochastic quantity, and cannot 
be used to label the Langevin equation, otherwise one would describe a physically different 
process. The time label must therefore be carefully specified. The second one is that, since H 
is related to the curvature of space-time, its stochasticity has to do with the one of space-time 
itself. We are thus a priori describing effective quantum gravitational effects, corresponding 
to the gravitational- and self-interactions of the inflaton field. The corresponding corrections 
should therefore remain small as long as the energy density of the inflaton field is small 
compared to the Planck scale. For this reason, it is convenient to define the dimensionless 
potential 


V _ H 2 

247r 2 Mp 1 ~ 87r 2 Mp! 


( 2 . 1 ) 


which we will make use of extensively in the following. Before turning to the calculation of 
the correlation functions, in this section, we show on a simple example why different time 
labels in the Langevin equation typically yield results that differ by oc v corrections. 

First, let us recast the stochastic process (1.1) through a Fokker-Planck equation, which 
governs the time evolution of the probability density P(<j>, N ) that ip = <f> at time N. In the 
ltd interpretation 4 5 [17, 35, 36], it reads 6 


d 

dN 


P (0) N) 


d 

d(j) 


V ' 
3 H 2 


P{<t>,N) 



' H 2 
8vr 2 


P(</>,N) 


(2.3) 


4 Hereafter, by “stochastic quantities”, we simply refer to realization dependent quantities, as opposed to 
quantities that are fixed for all realizations. 

5 More generally, the last term in Eq. (2.3) can be written in the form 


d_ I H 2a d 
dcj) \ 8tt 2 d(j> 


[n 2( i““ ) P(0,i V )] | 


( 2 . 2 ) 


with 0 < a < 1, where a = 0 corresponds to the Ito interpretation and a = 1/2 to the Stratonovich one [34]. 
However, analysis shows that keeping terms explicitly depending on a exceeds the accuracy of the stochastic 
approach in its leading approximation (1.1). In particular, corrections to the noise term due to self-interactions 
of small-scale fluctuations (if they exist) are at least of the same order or even larger. 

®Note also that we never use the “volume weighted” variant of Eq. (2.3) proposed as an alternative in 
Ref. [24] since then the resulting distribution is not normalizable: its integral over d(f> is time- or IV-dependent. 
Thus, it leads to probability non-conservation. Neither is it justified from the physical point of view, since it 
is based on the assumption that all Hubble physical volumes (“observers”) emerging from the expansion of a 
previous inflationary patch are clones of each other, while they are strongly correlated. 
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Now let us compare this equation with the one that would have been obtained if the Langevin 
equation had been written in terms of cosmic time t. Performing the simple change of time 
variable dlV = Hdt in Eq. (1.1), this is given by 


dp V' H 3 / 2 _ . , 

c it ~ ~3H + 2tt ^ ' 


(2.4) 


Here we use the notation p to stress the fact that, a priori, p does not describe the same 
stochastic process as p. The Fokker-Planck equation corresponding to Eq. (2.4) is given by 




TmAO 


+ 


dcj) 2 


gmAO 


(2.5) 


If H is taken to be a function of time only, independent of p, the H factors can be taken out 
of the derivatives with respect to cj> in Eqs. (2.3) and (2.5). In this case, it is straightforward 
to see that these two are perfectly equivalent through the change of time variable d N = Hdt, 
and that they describe the same stochastic process. On the contrary, if H explicitly depends 
on p, this is obviously no longer the case and P ^ P. 

This can be better illustrated by calculating the stationary distributions associated with 
these processes. Let P s t a t(</>) denote a stationary probability distribution for the stochastic 
process (1.1), or equivalently, (2.3). By definition, dP s ta.t{<j))/dN = 0, hence 


d_ 

dcj> 


V' 

2 -fstat (</*) + 


d_ 

dcj) 


\El P 

gyj-2 




( 2 . 6 ) 


which defines the probability current J. This current thus needs to be independent of cj> for 
a stationary distribution. In most interesting situations, it is actually 0. This is notably 
the case when the allowed values for cf> are unbounded. For example, if V(<f>) is defined up 
to i j) = oo, the normalization condition f P stat (cj))d(j) = 1 requires that P s ta,t(<j>) decreases 
at infinity strictly faster than |</>| _1 . In this case, both P s tat (0) and dP s ta,t(<j>) / dcj> vanish at 
infinity. From Eq. (2.6), J vanishes at infinity also, hence everywhere. This yields a simple 
differential equation to solve for -P s tat (</>)) and one obtains 


Tstat (4 > j OC 


1 


v( 4 >) 


exp 


1 


v{(j>) 


(2.7) 


Here, an overall integration constant, which makes the distribution normalized, J P s t a t(0)d</> = 
1, is omitted. Similarly, Eq. (2.5) can be written as dP/dt = dJ/dcj), and requiring that the 
current J vanishes gives rise to a differential equation for the stationary distribution P s ta.t((j>), 
which can easily be solved. One obtains 


P 


stat 



r 1 i 

3/2 

r 1 i 

oc 

y(4>). 

exp 



( 2 . 8 ) 


The two distributions are close, and the effects coming from the H(p) dependence are small, 
only in the regions of the potential where v <C 1. 

At this point, we are left with the issue of identifying the right time variable to work 
with. Actually, one can explicitly show [26, 27, 37] that N is the correct answer, and that 
it is the only time variable that allows the stochastic formalism to reproduce a number of 
results from QFT on curved space-times. We leave this discussion to appendix A, where we 
elaborate on existing results and show why, since we deal with metric perturbations, we must 
work with N. 
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3 Method 


Let us now review how correlation functions of curvature fluctuations can be calculated in 
stochastic inflation, and see which approach is best suited to the issue we are interested in. 

The problem can first be treated at linear order [38-40] by expanding the coarse grained 
field <p about its classical counterpart 4> c \, <p = (f> c \ + Here, recall that (f> c \ is the solution 
of the Langevin equation (1.1) without the noise term. The quadratic moment of can 
be calculated as in appendix A.2, see Eq. (A.23). It corresponds to the integrated power 
spectrum of the field fluctuations on super-Hubble scales, and can therefore be related [40] 
to the power spectrum V £ of curvature perturbations thanks to the relation 


V\ ~ 



(3.1) 


In this expression, the right hand side needs to be evaluated when the scale associated with 
the wavenumber k (at which the power spectrum is calculated) exits the Hubble radius. If 
one plugs the expression (A. 23) obtained in appendix A for using N as the time 

variable into Eq. (3.1), one obtains 


Vl ^ H2 ^ ci) 

C 87r 2 A/| 1 ei (^ c i) 


(3.2) 


As before, 4> c \ needs to be evaluated when the scale associated with the wavenumber at which 
the power spectrum is calculated exits the Hubble radius. The quantity e\ = — dH/(H 2 dt) 
is the first slow-roll parameter. At leading order in slow roll, it verifies ei = Mp l /2(v > /v) 2 . 
The above expression exactly matches the standard result [41, 42]. In order to get the first 
corrections to this standard result, one thus needs to go to higher orders in 5(f). Actually, one 
can show that no contributions arise at next-to-leading order, and that one needs to go at 
least to next-to-next-to-leading order. This renders the calculation technically difficult. This 
is why we will prefer to make use of non-perturbative techniques. In passing, let us stress 
that in Ref. [40], the Langevin equation is written and solved with t, whereas, as already said, 
the number of e-folds N must be used instead. This has important consequences. Indeed, if 
one makes use of cosmic time t and plugs the associated expression (A.29) for the quadratic 
moment of into Eq. (3.1), one obtains 


'jy ~ 

7 C - 


87r 2 Mpiei (0 c i) 


1 + 2 


H' (0d) l 2 [** 


' H(<j> ci) l 3 
H' (</> cl )_ 



(3.3) 


Here, we have adopted the same notation as in section 2 where a tilde recalls that not the 
same quantity is worked out and ( is not (. This result matches Eq. (2.11) of Ref. [40]. 
However, in this work, it is concluded that, because of the second term in the braces of 
Eq. (3.3), which is always negative, the amplitude of the power spectrum in the stochastic 
approach is in general reduced with respect to the standard result. One can see that such a 
statement is incorrect, since the extra term in Eq. (3.3) is simply due to not working with 
the correct time variable. This is why, if such an approach were to be followed and extended 
to higher orders, it would again be crucial to work with N as the time variable. 

Another strategy is followed in Refs. [43-45], where methods of statistical physics, such 
as replica field theory, are employed in a stochastic inflationary context. However, only the 
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case of a free test field evolving in a de Sitter or power-law background is investigated, while 
we need to go beyond the fixed background assumption in order to study the effects of the 
explicit H((p) dependence. This is why we cannot directly make use of this computational 
scheme in the present work. 

Finally, in Refs. [33, 46, 47], the 5N formalism is used to relate the curvature perturba¬ 
tions to the number of e-folds statistics. This is this last route that we chose to follow here, 
since it does not rely on any perturbative expansion scheme, and since it does not prevent 
us from implementing the explicit H(ip) dependence. In Ref. [47], numerical solutions are 
obtained for quadratic and hybrid potentials. In the present work, we derive fully analyti¬ 
cal and non-perturbative results that apply to any single-field potential, and which do not 
require a numerical solution of the Langevin equation. As a by-product, this allows us to 
prove, for the first time, that the standard results are always recovered in the classical limit, 
for any potential. 

3.1 The SN Formalism 

The 5N formalism [9, 48-52] is very well suited to addressing the calculation of correlation 
functions in stochastic inflation, since it relates the statistical properties of curvature pertur¬ 
bations to the distribution of the number of e-folds among a family of homogeneous universes. 
Let us first recall where this correspondence comes from and, as an example, how the scalar 
power spectrum is usually calculated in the associated formalism. 

Starting from the unperturbed flat Friedmann-Lemaitre-Robertson-Walker (FLRW) line 
element, ds 2 = — dt 2 + a 2 (t)5ij dxMad, deviations from homogeneity and isotropy can be 
included in a more general, perturbed metric, which contains some gauge redundancy. A 
specific gauge choice consists in requiring that fixed t slices of space-time have uniform 
energy density, and that fixed x worldlines be comoving. When doing so, and including 
scalar perturbations only, the perturbed metric in this gauge (which coincides in the super- 
Hubble regime with the synchronous gauge supplemented by some additional conditions 
fixing it uniquely) becomes [9, 53, 54] ds 2 = —dt 2 + a 2 (t)e 2 ^ x ^6ijdx‘ l dx : >, up to small terms 
proportional to gradients of Here, £ is the adiabatic (curvature) perturbation, which is 
time-independent in single-field inflation once the decaying mode can be neglected. The 
omission of tensor perturbations is justified by the fact that their amplitude is suppressed 
compared to the scalar ones by the small slow-roll parameter e\. This allows us to define a 
local scale factor a(t,x ) = a(t)e^ x \ Starting from an initial flat slice of space-time at time 
tin, the amount of expansion N(t,x ) = In [a(t, cc)/a(ti n )] to a final slice of uniform energy 
density is then related to the curvature perturbation through 

C{x) = N (t,x) - N 0 (t) = SN, (3.4) 

where No(t) = In [a(t)/a(tj n )] is the unperturbed amount of expansion. From this, an im¬ 
portant simplification arises on large scales where anisotropy and spatial gradients can be 
neglected, and the local density, expansion rate, etc., obey the same evolution equations as a 
homogeneous FLRW universe. Thus we can use the homogeneous FLRW solutions to describe 
the local evolution, which is known as the “quasi-isotropic” [55-58] or “separate universe” 
approach [51, 59, 60]. It implies that N(t,x) is the amount of expansion in unperturbed, 
homogeneous universes, so that £ can be calculated from the knowledge of the evolution of a 
family of such universes. Written in terms of the inflaton field <f)(x) = (f> + 5cf>(x), consisting 
of an unperturbed, homogeneous piece <f> and of a perturbation 5<f> originating from quantum 
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fluctuations, Eq. (3.4) gives rise to 


C (x) = N [p (t), cj) (as)] -N[p (t), <j>\ . 


(3.5) 


Here, N is to be evaluated in unperturbed universes from an initial epoch when the inflaton 
field has an assigned value 0 to a final epoch when the energy density has an assigned value 
p. Since the observed curvature perturbations are almost Gaussian, at leading order in 
perturbation theory, one has 

dN 

( (x) = 5N ~ -rr-5(j ). (3-6) 

t](p 

Here, N ((/>) is usually evaluated with the slow-roll, classical formula 


N(</>) 


1 

Mpi 


d(j> 


(3.7) 


Once ( is decomposed into Fourier components, (k = (27r)~ 3 / 2 f d 3 x((x) exp(ik ■ x), the 
power spectrum V( is defined from the quantum expectation value (CfcCfc') = (27t) 5 /(2/c 3 ) 
V^(k)6(k + k'). It can be expressed in terms of the power spectrum of 5(j) (defined by similar 
relations) thanks to Eq. (3.6). For quasi-de Sitter inflation, and when the curvature of the 
inflaton potential is much smaller than H, on super-Hubble scales, the latter is given by [61] 
Vs(j>{k) — H 2 (k)/4:n 2 , where H(k ) means H evaluated at the time when the k mode crosses 
the Hubble radius, i.e. when aH = k. Together with Eq. (3.7), one therefore obtains 


n = 


'H(k) l 2 1 


2vr 


2Mp,ei (k) ’ 


(3.8) 


which is the same as Eq. (3.2) and which matches the standard result [41, 42], 

A fundamental remark is that in the above usual calculation, the quasi-isotropic (sep¬ 
arate universe) approximation is assorted with the assumption that on super-Hubble scales, 
the evolution of the inflaton field is governed by its classical equation of motion (3.7). The 
stochastic dispersion in the number of e-folds thus only comes from the field dispersion at 
Hubble crossing 5(f >*. In most cases, this is a good approximation for the following reason. 
From the Langevin equation (1.1), one can see that during the typical time scale of one e-fold, 
the classical drift of the inflaton field is of the order A</> c i = V'/(3H 2 ) = \f2e\M Pi . while the 
quantum kick is of the order Ac/> qu = H/(2n). This allows us to define a rough “classicality” 
criterion A<p qu /A(p c i that assesses the amplitude of the stochastic corrections to the classical 
trajectory. Making use of Eqs. (3.8), this ratio can be expressed as 


A^qu 

A^ci 


\/A. 


(3.9) 


which is valid for single-field slow-roll models of inflation with canonical kinetic terms. Since 
V ~ 2 x 10~ 9 for the modes observed in the Cosmic Microwave Background (CMB), stochastic 
effects are already small when these modes cross the Hubble radius. If one further assumes 
that ei monotonously grows toward 1 during the last stages of inflation, V( oc H 2 /e\ decreases 
(since H can only decrease) and one is therefore ensured that the stochastic corrections to 
the inflaton trajectory remain small. 

However, they are two caveats to this line of reasoning. The first one is that, as we will 
show below, A(j> qu /A(j> c \ is not the correct way to assess the importance of stochastic effects 
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and one should use instead another classicality criterion that we will derive. The second one 
is that, in some situations, becomes tiny or even vanishes in some transient phase between 
the Hubble exit time of the observed modes and the end of inflation. This is the case, for 
example, when the potential has a flat inflection point, such as in MSSM inflation [62-64] or 
as in punctuated inflation [65, 66]. Another situation of interest is when inflation does not 
have a graceful exit but ends due to tachyonic instability involving an auxiliary field, like in 
hybrid inflation [67, 68], or by brane annihilation in string-theoretical setups [69, 70]. In such 
cases, €i can decrease and the last e-folds of inflation may be dominated by the quantum 
noise. It is therefore important to study the dispersion 5N arising not only from 5<j)* but 
from the complete subsequent stochastic history of the coarse grained field. 

Note also that in these expressions, £ need not be small as was shown in Refs. [9, 51, 71] 
[note, however, that ( is defined up to a constant due to an arbitrary possible rescaling 
of a(f)], thus, 5N need not be small, too. As follows from the quasi-isotropic (separate 
universe) approach, the condition for inflation to proceed is only that H <C M P1 . On the 
other hand, if P^(fc) ~ H/(M Ply /e i) exceeds unity (the so called regime of “eternal inflation”), 
then the Universe loses its local homogeneity and isotropy after the end of inflation, but 
not immediately. This occurs much later than the comoving scale a(t)/k at which this 
inhomogeneity occurs crosses the Hubble radius H^ 1 second time. Thus, in the scope of 
the inflationary scenario may well exceed unity at scales much exceeding the present 
Hubble radius. The stochastic inflation approach provides us with a possibility to obtain 
quantitatively correct results in this non-linear regime, too. 


3.2 Computational Programme 

This is why we now generalize this approach to a fully stochastic framework. For a given 
wavenumber k, let (j>*(k) be the mean value of the coarse grained field when k crosses the 
Hubble radius. If inflation terminates at </> en d, let Af(k) denote the number of e-folds realized 
between </>*(&) and 0 en d- Obviously, Af is a stochastic quantity, and we can define its variance 

5Af 2 (k) = (Af 2 (k)) - (Af ( k )) 2 . (3.10) 


It is related with the curvature perturbation 5N of Eq. (3.6) in the following manner. Since 
5AT is computed between two fixed points <^>* (A;) and cWd, it receives an integrated contri¬ 
bution from all the modes crossing the Hubble radius between these two points, and one 
has 


5N 2 (k) 



'PSN (k) 


d k 
k 


/*ln k end 

An k end - (Af) (1—ei*+---) 


VsNdN. 


(3.11) 


Here we have used the relation (Af(k)) = ln(a en d/a*) = ln(fc Rn d/fc)(l + ei* + ■■•), where 
ei* + ■ • • stand for slow-roll corrections that we do not need to take into account at leading 
order in slow roll. One then has 


Pc (*) = Pmc (k) 


ddAr 2 

d W) 


(AC)=ln(fe end /fe) 


(3.12) 


In the same manner, the third moment of the number of e-folds distribution, 


6Af 3 (k) = ((Af-(Af)) 3 ), 


(3.13) 


receives a double integrated contribution from the local bispectrum and one has B £ oc 
d 2 5Af 3 /d{Af} 2 . The local /nl parameter, measuring the ratio between the bispectrum and 
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Figure 1. Sketch of the dynamics solved in section 3. The inflaton is initially located at (/* and 
evolves along the potential V(<j>) under the stochastic Langevin equation (1.1), until it reaches one of 
the two ending values </>i or 0 2 - The left panel is an example where inflation always terminates by 
slow-roll violation, while the right panel stands for a situation where one of the ending points, <f> 2 , 
corresponds to where V ~ Mp, above which inhomogeneities prevent inflation from occurring. 


the power spectrum squared, is then given by 

(3.14) 

where 5/72 is a conventional historical factor. Analogously, the trispectrum is related to the 
third derivative of 5AT 4 with respect to (J\f), and so on and so forth. 

The computational programme we must follow is now clear. For a given mode k, we 
first calculate (j>*(k) (this sets the location of the observational window). We then consider 
stochastic realizations of Eq. (1.1) that satisfy ip = (j>*(k) at some initial time,' and denote by 
A f the number of e-folds that is realized before reaching </> en d- Among these realizations, we 
calculate the first moments of this stochastic quantity, (AT), (Af 2 ), (AT 3 ), etc. We finally apply 
relations such as Eqs. (3.12) and (3.14) to obtain the power spectrum, the non-Gaussianity 
local parameter, or any higher order correlation function. 

3.3 First Passage Time Analysis 

In what follows, this calculation is performed using the techniques developed in “first passage 
time analysis” [72, 73], which was applied to stochastic inflation in Ref. [17]. We consider the 
situation sketched in Fig. 1, where the inflaton is initially located at <^* and evolves in some 
potential V (4>) under Eq. (1.1). Because any part of the potential can a priori be explored, 

7 This calculation therefore relies on a specific choice of initial (in fact, pre-inflationary) conditions, since 
all trajectories emerge from </>» at initial time. In principle, other choices could be made, even if most physical 
quantities (in particular, perturbations during the last 60 e-folds) do not depend on them. 


5 d 2 5Af 3 fd5AT 2 y 2 
fNL ~ 72d(A02 Vd(A0/ ’ 














here we consider two possible ending points, 4>i and <f> 2 , located on each side of 1 j>*. If the 
potential is, say, of the hilltop type (left panel), (f> 1 and (j )2 can be taken at the two values 
where inflation has a graceful exit, on each side of the maximum of the potential. If, on the 
other hand, a flat potential extends up to 0 = 00 (right panel), one of these points, say <f> 2 , 
can be taken where V becomes super-Planckian and inhomogeneities prevent inflation from 
occurring. In such cases, the precise value of </>2 plays a negligible role, as we will show in 
section 3.3.1. Let A f be the number of e-folds realized during this process. 

Before proceeding with the calculation of the J\f moments, a first useful result to establish 
is the Ito lemma, which is a relation verified by any smooth function / of p. The Taylor 
expansion of such a function at second order is given by f {p + dp) = / {p) + f (p) dp + 
f" (p) /2dp 2 + O (dy? 3 ). Now, if p is a realization of the stochastic process under study, dp 
is given by Eq. (1.1) and at first order in dIV, one obtains 

df [p ( N )] = /' \p (IV)] y/2v [p(N)]M P £ (N) dN 

~f [</? {N)] 1>[p(N)\ M * ldN + M * J " [</? {N)] V ^ {N)] dN ' (3 - 15) 

Integrating this relation between IV = 0 where p = </>* and N = AT where p = <p\ or </> 2 , one 
gets the Ito lemma [74] 

rJ\f 

f (0i or 4 > 2 ) - / (0*) = / f'[p (IV)] sj2v [p (IV)]M P1 £ ( N) dN 

Jo 

+J* {<r [p m v [p (iv)] - / [0 (iv)] div, 

(3.16) 

which we will repeatedly make use of in the following. 


3.3.1 Ending Point Probability 

As a first warm-up, let us calculate the probability pi that the inflaton field first reaches the 
ending point located at 4>i [i.e. (f> (, M ) = (j) 1 ], or, equivalently the probability P 2 = 1 — Pi that 
the inflaton field first reaches the ending point located at (j )2 [i.e. cj) (J\f ) = 4> 2 ] ■ This will also 
allow us to determine when the ending point located at 4 >2 plays a negligible role. 

First of all, let i/j (p) be a function of the coarse grained field that can be expressed as 


'!>(&) 


h(p) -hjfc) 
h(<f>i) -h((j)2) ’ 


(3.17) 


where h (p) will be specified later. By construction, one has tj) (^>i) = 1 and ijj (^> 2 ) = 0- This 
implies that the mean value of evaluated at p(N) is given by {ip [p (J\T)]) = piV’(^i) + 
P 2 1/) (</> 2 ) = Pi- The idea is then to find an h (p) function that makes easy the evaluation of 
the left hand side of the previous relation, so that we can deduce p \. In order to do so, let 
us apply the Ito lemma (3.16) to h (p). If one requires that the integral of the second line of 
Eq. (3.16) vanishes, that is, 

h" ( P) v (p) = ti (p) ’ ( 3 - 18 ) 

v{p) 

one obtains 


rJ\f 

h [p (AT)] - h (<f>i n ) = / ti [p (IV)] sj2v [p (IV)]M P1 £ (IV) dA^. (3.19) 

Jo 
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Because h and ip are linearly related, see Eq. (3.17), the same equation is satisfied by ip. When 
averaged over all realizations, * * 8 its right hand side vanishes. One then obtains (ip [ip (JV)]) = 
V 7 (0*), which is the probability p\ one is seeking for. All one needs to do is therefore to 
solve Eq. (3.18) to obtain h(ip), to plug the obtained expression into Eq. (3.17) to derive 
ip (fa), and finally to evaluate this function at fa. A formal solution to Eq. (3.18) is given by 
h (fa) = A exp [— 1/v (x)] dx, where A and B are two integration constants that play no 
role, since they cancel out when calculating ip thanks to Eq. (3.17). Indeed, the latter gives 
rise to 



(3.20) 


and a symmetric expression for p 2 . 9 

A few remarks are in order about this result. First, one can check that, since fa lies 
between fa and fa, the probability (3.20) is ensured to be between 0 and 1. Second, one 
can also verify that when fa = fa, p\ = 1, and when fa = fa, p\ = 0, as one would expect. 
Third, in the case depicted in the right panel of Fig. 1, in the limit where fa —> oo, one is sure 
to first reach the ending point located at fa, that is, p 2 = ff* e ~ 1 / v / fj’ 2 e~ l / v = 0. Indeed, 
the numerator of the expression for p 2 is finite, since a bounded function is integrated over a 
bounded interval. If the potential is maximal at fa, and if it is monotonous over an interval 
of the type [fa, fa[, its denominator is on the contrary larger than the integral of a function 
bounded from below by a strictly positive number, over an unbounded interval [fa, (poo[- This 
is why it diverges, and why p 2 vanishes. This means that if fa is sufficiently large, its precise 
value plays no role, since inflation always terminates at fa = fa n d- 


3.3.2 Mean Number of e-folds 

Let us now turn to the calculation of the mean number of e-folds (A f). As above, we want 
to make use of the Ito lemma (3.16). To do so, let us define f(fa) as the solution of the 
differential equation 

= (3.21) 

with boundary conditions / (fa) = / (fa) = 0. Such a solution will be explicitly calculated 
in due time. For now, it is interesting to notice that when this is plugged into the Ito 
lemma (3.16), the first term of the left hand side, f (fa or fa), vanishes, and the second 
integrand of the right hand side is —1. Thus, the Ito equation can be rewritten as 


Af 


f (fa) + 



f Yp (N)] V 2v [v(N)\Mr£ (N ) d N . 


(3.22) 


s The fact that the averaged integral in the right hand side of Eq. (3.19) vanishes is non-trivial since both 

the integrand and the upper bound are stochastic quantities, but this can be shown rigorously (see e.g. p. 12 

of Ref. [73]). 

9 This is in agreement with Eq. (29) of Ref. [17], derived in the case where H is constant, hence 1 E 1 « 
vP 1 — (v — i>*) vZ 2 , where fa and fa lie at ±00 correspondingly, and where the initial condition for Eq. (2.3) 
is chosen to be P(<j), 0) = S{4> — fa). 
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Figure 2. Integration domain of Eq. (3.24) when evaluated at ip = (f> 2 , in the case 4>\ < </>2 (the 
opposite case proceeds the same way). The discrete parameter x is integrated between </> 1 and <fi 2 , while 
y varies between x and (f>. The resulting integration domain is displayed in green. When x < <p, one 
has dxdy > 0 and one integrates a positive contribution to the mean number of e-folds. Conversely, 
when x > (f>, one has Axdy < 0 and one integrates a negative contribution. This is necessary in order 
for the overall integral to vanish. This is why cf> must lie between cf> 1 and (f> 2 . 


By averaging over realizations, one obtains 10 


(A0 = /(</>*) • 


(3.23) 


What one needs to do is therefore to solve the deterministic differential equation (3.21) with 
the associated boundary conditions, and to evaluate the solution at One obtains 


f(<p) 



d y 1 

Mpi v ( y ) eXP 



1 

v ( x ) 


(3.24) 


where cj) is an integration constant set to satisfy the condition /(</> 2 ) = 0. There is no generic 
expression for it, 11 but one can be more specific. First of all, as can be seen in Fig. 2, ^ must 
be such that, when / is evaluated at cj) 2 , the integration domain of Eq. (3.24) possesses a 
positive part and a negative part, which are able to compensate for each other. This implies 
that (f> must lie between cf> 1 and </> 2 . A second generic condition comes from splitting the 
x-integral in Eq. (3.24) into dx = dx + dx. The first integral vanishes because 

10 Here again, since both the integrand and the upper bound are stochastic quantities, it is non-trivial that 
the integral in the right hand side of Eq. (3.22) vanishes when averaged, but it can be shown rigorously. 

11 Alternatively, one can write Eq. (3.24) in the explicit form [17] 


to 

a 

*5 

dx 

1 

1 

1 

Tfpi 

Jy Af Pl 

/ \ 6xp 

v(y) 

y(y) 

v(x) 


\6(x - x,) - pi] , 


where pi is given by Eq. (3.20) and, in the configuration of Fig. 1, 9(x — x*) = 1 when x > x* and 0 otherwise. 
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/( 02 ) = 0 , which means that in order for / to be symmetrical in 0 i <->• 02 , 0 ( 0 i, 02 ) must 
satisfy this symmetry too, that is to say, 0 (01, 02) = 0 (02, 0i). Third, in the case where the 
potential is symmetric about a local maximum 0 max close to which inflation proceeds, the 
integrand in Eq. (3.24) is symmetric with respect to the first bisector in Fig. 2. The two green 
triangles must therefore have the same surface, which readily leads to 0 = 0 max • Fourth, 
finally, in the case displayed in the right panel of Fig. 1, if 0 2 is sufficiently large, we have 
established in section 3.3.1 that P 2 = 0 and the quantity we compute is the mean number 
of e-folds between 0* and 0i = 0 en d- For explicitness, let us assume that v' > 0 (the same 
line of arguments applies in the case v' < 0). Inflation proceeds at 0 < 02- In the domain of 
negative contribution in Fig. 2, the argument of the exponential in Eq. (3.24) is positive. As 
a consequence, if 0 is finite and 02 —»• 00 , the negative contribution to the integral is infinite 
while the positive one remains finite, which is impossible. In order to avoid this, one must 
then have 0 = 0 2 . In practice, almost all cases boil down to one of the two previous ones 
and 0 is specified accordingly. Combining Eqs. (3.23) and (3.24), one finally has 12 



d y 1 

Mpi v (■ y ) 6XP 



1 

v (x) 


(3.25) 


This quantity is plotted for large and small field potentials in Fig. 3, where it is compared with 
the results of a numerical integration of the Langevin equation (1.1) for a large number of 
realizations over which the mean value of M is computed. One can check that the agreement 
is excellent. 


Classical Limit 

Let us now verify that the above formula boils down to the classical result (3.7) in some 
“classical limit”. This can be done by performing a saddle-point approximation of the 
integrals appearing in Eq. (3.25). Let us first work out the y-integral, that is to say, 
J x> dy/v{y) exp[l/u(y)]. Since the integrand varies exponentially with the potential, the strat¬ 
egy is to evaluate it close to its maximum, i.e. where the potential is minimum. The potential 
being maximal at 0 in most cases (see the discussion above), the integrand is clearly maxi¬ 
mal 13 at x. Taylor expanding \/v at first order around x. 1 /v(y) — l/v(x)—v'(x)/v 2 (x)(y—x), 
one obtains, after integrating by parts, 14 dy/v(y) exp [l/v(y)\ — v{x)/v'(x) exp [\/v{x)\. 
Plugging back this expression into Eq. (3.25), one finally obtains 


(AO L = 


f A* dx v(x) 

0end M P1 V '( X ) ’ 


(3.26) 


which exactly matches the classical result (3.7). The classical trajectory thus appears as a 
saddle-point limit of the mean stochastic trajectory, analogously to what happens e.g. in the 
context of path integral calculations. 

This calculation also allows us to identify under which conditions the classical limit is 
recovered. A priori, the Taylor expansion of 1/v can be trusted as long as the difference 
between l/v(x) and 1 /v(y) is not too large, say |1 /v(y) — l/v{x)\ < R, where R is some 

12 This is again in agreement with Eq. (35) of Ref. [17] if H is constant and 0* = 0 = 0, while 0 en d = 00 . 

13 Strictly speaking, this is only true if the potential is a monotonous function of the field, but this is most 
often the case in the part of the potential that is relevant to the inflationary phase. 

14 Since v(0) v(x) and if v is monotonous, one can also show that exp [— v'{x)/v 2 (x){4> — a:)] is exponen¬ 
tially vanishing and this term can be neglected. 


13 
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Figure 3. Mean number of e-folds (AT) realized in the large field V oc </> 2 (left panel) and small field 
V oc 1 — (j) 2 /y? (where ji = 20M P1 , right panel) potentials, as a function of </>». The location </>®° refers 
to the value of (/)* for which the classical number of e-folds N c \ = 50 and <j > en( j is where ei = 1. In both 
panels, the overall mass scale in the potential is set to the value that fits the observed amplitude of 
the power spectrum ~ 2 x 10~ 9 , when the latter is calculated with the classical formula (3.8), 50 
e-folds before the end of inflation. The green line corresponds to the analytical exact result (3.25), 
and the red circles are provided by a numerical integration of the Langevin equation (1.1) for a 
large number of realizations over which the mean value of M is computed. The orange dashed 
line corresponds to the classical limit (3.26). The top axes display v and the classicality criterion 
\2v — v"v 2 /v' 2 \. The yellow shaded area stands for v > 1, where inhomogeneities are expected to 
prevent inflation from occurring and our calculation cannot be trusted anymore. 


small number. If one uses the Taylor expansion of 1/v at first order, this means that | y — 
x\ < Rv 2 /v'. Requiring that the second order term of the Taylor expansion is small at the 
boundary of this domain yields the condition \2v — v"v 2 /v' 2 \ 1. For this reason, we define 

the classicality criterion 


Vcl = 


2v - 


v"v 2 


(3.27) 


This quantity is displayed in the top axes in Fig. 3 and one can check that indeed, the 
classical trajectory is a good approximation to the mean stochastic one if and only if 7] c \ <C 1. 
In the following, we will see that r] c \ is the relevant quantity to discuss the strength of the 
stochastic effects in general and in section 4.3, we will further discuss the physical implications 
of Eq. (3.27). 

For now, and for future use, let us give the first correction to the classical trajectory. 
This can be obtained going one order higher in the saddle-point approximation, that is to 
say, using a Taylor expansion of 1/v at second order. One obtains 


m 


7?cl<l 


rdx v(x) 

4end M Pl V '( X ) . 


v” (x) v 2 (x) 

l + v(x) -+ 


v' 2, (i X ) 


(3.28) 
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where the dots stand for higher order terms. In the brackets of Eq. (3.28), the two last terms 
stand for the first stochastic correction and one should not be surprised that, in general, 
when r/d is small, it is small. It is also interesting to notice that it is directly proportional 
to dei/dN. When e\ increases as inflation proceeds, the stochastic leading correction is 
therefore positive and the stochastic effects tend to increase the realized number of e-folds, 
while when e± decreases as inflation proceeds, the correction is negative and the stochastic 
effects tend to decrease the number of e-folds, at least at linear order. 


3.3.3 Number of e-folds Variance 

Let us now move on with the calculation of the dispersion in the number of e-folds, defined 
in Eq. (3.10). If one squares Eq. (3.22), and takes the stochastic average of it, one obtains 15 

(M 2 ) = f 2 (fa) + 2 Ml f 2 [<f> (N )] v [0 (IV)] d N^j . (3.29) 

In order to make use of the Ito lemma, let then g(fa be the function defined by 

d" (<t>) v (4>) ~ a' (</>) V-yrl = -2 f 2 (fa v (fa , (3.30) 

v(fa 

where / is the function defined in Eq. (3.24). When the Ito lemma (3.16) is applied to g (fa, 
if one further sets g (fa) = g (fa) = 0, one obtains 

g (fa) = 2 Ml f [0 ( N)} v [</> (IV)] d 

= (AT 2 ) - f 2 (fa) = (M 2 ) - (AT ) 2 , (3.31) 


where the second equality is just a consequence of Eq. (3.29) and where the third equality is 
just a consequence of Eq. (3.23). Therefore, one just needs to solve Eq. (3.30) with boundary 
conditions g (fa) = g (fa) = 0 and to evaluate the resulting function at fa in order to obtain 
5Af 2 = g(fa). The differential equation (3.30) can formally be integrated, and one obtains 


a (fa) = 2 / dx 

J <f> 1 


(- 02 ( 01 , 02 ) 


dy/ /2 (y) exp 


_v(y) v(x)_ 


(3.32) 


where fa ( fa , fa) is an integration constant that must be chosen in order to have g (fa) = 0. 
One can show that it satisfies the four properties listed in section 3.3.2 for cf> and can therefore 
be specified in the same manner. With fa = <^ e nd, one then has 


/“0end 

5M 2 = 2 / dx 


dyf ' 2 ( y ) exp 


'02 


1 


1 


v(y) v(x) 


(3.33) 


Classical Limit 

As was done for the mean number of e-folds in section 3.3.2, let us derive the classical limit of 
Eq. (3.33). Obviously, in the classical setup the trajectories are not stochastic and 5J\f 2 = 0, 
and what we are interested in here is the non-vanishing leading order contribution to 5M 2 

15 This is again a non-trivial result since both the integrand and the upper bound of the integral appearing 
in Eq. (3.22) are stochastic quantities, but, as before, it can be shown in a rigorous way. 
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in the limit ?y c i < 1. As before, the y-integral can be worked out with a saddle-point ap¬ 
proximation, and one obtains 16 ff 2 dyf' 2 (y) exp [1 /v(y)\ — v 4 ( x) /v' 3 (;x ) exp [\/v (x)] /Mp r 
Plugging back this expression into Eq. (3.33), one obtains 


5N 2 


cl 


2 r+* v A {x) 

M Pi4end X V ,3 W 


(3.34) 


Finally, and for future use again, let us give the first correction to this classical limit. Going 
one order higher in the saddle-point approximation, one obtains 


5N 2 


T?cl<l 


2 ft* v A (x) 

M Pi4end X V' 3 (X ) 


1 + 6v (. x ) 


v 2 (x) v" ( x ) 
v' 2 (x) 



(3.35) 


3.3.4 Number of e-folds Skewness and Higher Moments 

In the same manner, if one denotes the third moment of the distribution of number of e- 
folds by m ( q !>*) = 5AT 3 defined in Eq. (3.13), one can show that m((j)) is the solution of the 
differential equation m" — m'v'/v 2 = —6fg' that obeys m(<j) i) = rnifa) = 0. As before, 
taking cf> \ = ^) ent i and = <f>oc, one obtains 


5AT 3 = 6 



dyf'(y)g'(y)ex p 



1 

V (x) 


(3.36) 


where (f>$ can be set as </>. Similarly to above, making use of Eqs. (3.28) and (3.35), a 
saddle-point approximation of this integral leads to the classical limit 


5N 3 


»7cl<l 


12 [** , V 7 


^1 + 14u - 11 


v 2 v" 



(3.37) 


Let us finally explain how the same procedure can be iterated and higher order moments 
can be calculated. Let us denote the p th momentum of the number of e-folds distribution by 


= 5Af p = ((J\f - (A f)) p ) , (3.38) 

where, by convention, we set <7o = 1 and <ti = 0. As above, one can recursively show that a p 
is the solution of the differential equation 

a p ~ a p^2 = ~ 2 Pf' a 'p -1 -PiP- 1 ) ? 2 o v - 2 (3.39) 

satisfying a p i) = a p (fo) = 0. On then has 

/•‘/’end rx 

a P {<!>*) = dz / d y 

J (f)* J (f)p 

(3.40) 

When p = 2, this yields the variance (3.33); when p = 3, the skewness (3.36) is obtained; 
when p = 4, the kurtosis could be derived as well, and so on and so forth for any moment. 

16 In the r/d <S 1 limit, / is close to the classical trajectory (3.26) as shown in section 3.3.2, and one can 
take /' (y) ~ v (y) /V (y) M~ 2 . 


2 pf (y) cr'p -1 iy) + P (p - l) f 2 iy) (Jp -2 (y) 


exp 


1 

v(y) 


l 

V (x) 
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Figure 4. Scalar power spectrum Vq for the large field V oc (j) 2 (left panel) and small field V oc. 
1 — (f> 2 / fj, 2 (where fi = 20M P1 , right panel) potentials, as a function of </>*. The conventions are the 
same as in Fig. 3. The green line corresponds to the analytical exact result (4.1), and the orange 
dashed line to the classical limit (4.3). 


4 Results 

We are now in a position where we can combine the intermediary results of the previous 
sections to give explicit, non-perturbative and fully generic expressions for the first correla¬ 
tion functions of curvature perturbations in stochastic inflation. We first derive the relevant 
formulas and their classical limits, before commenting on their physical implications in sec¬ 
tion 4.3. 


4.1 Power Spectrum 

Following the programme we settled in section 3.2, if one plugs Eqs. (3.25) and (3.33) into 
Eq. (3.12), one obtains = </(0*)// / (</’*)> that is, 


V c 


= 2 


f [t dx 

1 


i 

1 

1 

[J*. M P1 

v (x] 

exp 

V ( X) 

V (</>*)_ 




dx J 

r<t>00 ( \y 

1 

i 

1 

\ 2 

> exp 

1 

1 

M P1 1 

Jx ^P1 

/ \ 

v{y) 

y(y) 

v (x) 

V (x) 

V (</>*)_ 


(4.1) 


In this expression, Vq (0*) stands for the power spectrum calculated at a scale k such that 
when it crosses the Hubble radius, the mean inflaton field value is <^>*. This formula provides, 
for the first time, a complete expression of the curvature perturbations power spectrum 
calculated in stochastic inflation. It is plotted for large and small held potentials in Fig. 4. 
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From this, a generic expression for the spectral index can also be given. Since, at leading 
order in slow roll, d/dln(k) — —dcj)/d(J\f) x d/dq i>, one has 


n s = 1 — 


<?"(</>) , f" {4>) 


+ 




(4.2) 


Here, for conciseness, we do not expand this expression in terms of integrals of the potential, 
but it is straightforward to do so with Eqs. (3.24) and (3.32). 

Classical Limit 

Before commenting further on the physical implications of the above result, let us make sure 
that in the classical limit, g c \ <C 1, the standard formula is recovered. Combining Eqs. (3.12), 
(3.26) and (3.34), one has 

2 ^ 


0 = 


Ml v' 2 


(<t>* 


(4.3) 


which exactly matches the usual result (3.8) at leading order in slow roll. This fully generalizes 
the work of Ref. [46] , where this agreement is shown but only in the case where the Hubble 
parameter varies linearly with <f> and if the noise has constant amplitude. Here we have 
extended this result to any potential is single-field slow-roll inflation, and included dependence 
of the noise amplitude on the coarse grained field. In the same manner, making use of 
Eqs. (3.26), (3.34) and (4.2) together, one obtains n s | cl = 1 — M 2 X 3 {v'/v) 2 — 2v"/v , which 
again matches the standard slow-roll result n s = 1 — 2ei — £2 where £2 = dlnei/dlV is the 
second slow-roll parameter, since at leading order in slow roll, one has £2 = 2M 2 x {v' 2 /v 2 — 
v"/v), and e\ has been given above. 

Let us now derive the leading order corrections to these “classical” results. This can be 
done making use of Eqs. (3.28) and (3.35) in Eqs. (4.1) and (4.2). For the power spectrum, 
one obtains 


ni««i (w = nic 

while for the spectral index one gets 


1 + 5v (<)>*) — 4 


v 2 ((j>*)v" (0*) 


„/2 


( ( p* 


(4.4) 


n s \ 


»7cl<l 


(</>* 


^slcl 1 


+ M 2 


U /2 (0*) v" z ((j)*) V (</>*) V (4>*) v'" (cj>*) 


„//2 


3v" ((/>*) — 2 


V (</>*) 


- 6 


„/2 


(0* 


+ 4- 


v' (cj>* 


(4.5) 


4.2 Non-Gaussianity and Higher Moments 


The local non-Gaussianity parameter can be calculated in the same manner, and Eq. (3.14) 
gives rise to 


5 fm" f"m'\ 
!nL ~72 ■ 


(4.6) 


For conciseness, this expression is not expanded in terms of integrals of the potential, but it 
is straightforward to do so with Eqs. (3.24), (3.32) and (3.36). 

Here also, we need to make sure that in the classical limit, the standard result is recov¬ 
ered. Combining Eqs. (4.6), (3.28), (3.35) and (3.37), one obtains 


5 

/nl| t , c1 «i = 


/2 // / /2 // /// // 2 ' 
V A V / V t V V 4 V 


6^ - 4- + v ( 25^ - 34- - 10^— + 24^ ] + O {v 2 ) 


(4.7) 
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The first two terms in the brackets match the usual result [75]. In contrast, it is important 
to stress that within the usual 5N formalism, the standard result cannot be obtained be¬ 
cause of the intrinsic non-Gaussianity of the fields at Hubble exit [75, 76]. Such effects are 
automatically taken into account in our formalism, which readily gives rise to the correct 
formula. 

Obviously, one can go on and calculate any higher order correlation function with 
Eq. (3.40). However, with the power spectrum and non-Gaussianity local parameter at 
hand, we already are in a position where we can draw important physical conclusions. 

4.3 Discussion 

A first important consequence of Eqs. (4.1) and (4.6) is the correctness of their classical 
limits. They show the validity of our computational programme for calculating correlation 
functions in general. This may be particularly useful for investigating other cases than single¬ 
field slow-roll inflation, especially when the standard procedure is difficult to follow. Indeed, 
our method can easily be numerically implemented, and it could then be applied to more 
complicated scenarios such as multi-field inflation where it has been shown [77] that the SN 
formalism retains reliable, modified kinetic terms where the stochastic inflation formalism has 
been generalized [78-81], etc. In particular, it is well suited to situations where stochastic 
effects dominate the inflationary dynamics in some parts of the potential [47, 68, 82] and 
where one must take the stochastic effects into account. 

Let us mention that within the CMB observable window, corrections to the classical 
results are always small, since one has 

Vd - T( (d + j) • (4.8) 

More precisely, Eqs. (4.4) and (4.7) can be recast as V$\ ^ ~ V $\ d [l + V(\ A (ei + £ 2 )] 

and /nl|jj c1 «i - /nl| c i - T 2 Vc\ d (38ef + fr 2 ei + §e 2 e 3 - f e|), where e 3 = dlne 2 /dA r is 
the third slow-roll parameter. For the scales of astrophysical interest today, in standard 
single-field slow-roll inflation, these corrections are therefore tiny. 

However, even if the stochastic effects within the CMB observable window need to be 
small, let us stress that the location of the observable window along the inflationary potential 
can be largely affected. This notably happens when the potential has a flat region between 
the location where the observed modes exit the Hubble radius and the end of inflation, as is 
the case e.g. in hybrid inflation or in potentials with flat inflection points. 

Another point to note is that, contrary to what one may have expected, the corrections 
we obtained are not controlled by the ratio A (f>q U /A.(j) c \ extensively used in the literature, but 
by the classicality criterion r) c \ derived in Eq. (3.27). This has two main consequences. 

First, rj c \ has dimension v , which means that it is Planck suppressed. 1 ' This makes 
sense, since, as noted in section 2, some of the corrections we obtained physically correspond 
to the self- and gravitational interactions of the inflaton field. 11. This is why it can be 

17 This remark also sheds some new light on the old debate [83-85] whether quantum gravitational correc¬ 
tions should affect inflationary predictions through powers of <f>/M P 1 or V/M p,. This analysis reveals V/M Pl 
corrections only, regardless of the value of <f>/M P 1 . 

18 For this reason, one may think that performing the calculation in Fourier space as we did does not allow us 
to properly account for self-interaction effects and that a real space calculation should be carried out instead. 
However, since the separate Universe approximation is exponentially well verified on large scales, this is not 
the case. Making use of the same formalism as in Ref. [25] , we have indeed explicitly checked that performing 
the calculation in real space leads to the same results as the ones presented here. 
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Table 1. Classicality criterion r] c \ defined in Eq. (3.27) for a few types of inflationary potentials. 
Except for “large field”, the expression given for r] c \ is valid close to the flat point of the potential. 


useful to compare our results with loop calculations performed in the literature by means 
of other techniques. In particular, the self-loop correction to the power spectrum is derived 
in Ref. [86], and graviton loop corrections are obtained in Ref. [87] (for a nice review, see 
also Ref. [88]). A diagrammatic approach based on the 5N formalism is also presented in 
Ref. [89] where the power spectrum and the bispectrum are calculated up to two loops. In all 
these cases, the obtained corrections are of the form V 1 Ioop = V^ Tee (l + aV^ ree e 2 N). Here, 
a is a numerical factor of order one that depends on the kind of loops one considers, and e 2 
stands for second order combinations of slow-roll parameters. When the number of e-folds N 
is of the order 1/e, this is exactly the kind of leading corrections we obtained. This feature 
is therefore somewhat generic. Obviously, it remains to understand which loops exactly our 
approach allows one to calculate, and how our results relate to the above mentioned ones. 
We leave it for future work. 

Second, rj c \ contains 1/v' 2 terms. This means that, even if v needs to be very small, 19 
if the potential is sufficiently flat, r] c \ may be large. In table 1, we have summarized the 
shape of ?y c i for different prototypical inflationary potentials. For large field potentials, rj c \ 
is directly proportional to v. This is why, in the left panels of Figs. 3 and 4, departure of 
the stochastic results from the standard formulas occur only when 1, in a regime where 
our calculation cannot be trusted anyway. However, for potentials with flat points, different 
results are obtained. If the flat point is of the hilltop type, r/ c \ diverges at the maximum of 
the potential. This is why, in the right panels of Figs. 3 and 4, even if v saturates to a small 
maximal value, the stochastic result differs from the classical one close to the maximum of the 
potential. However, this happens many e-folds before the scales probed in the CMB cross the 
Hubble radius, that is to say, at extremely large, non-observable scales. The same conclusion 
holds for plateau potentials (either of the polynomial or exponential type) where stochastic 
effects lead to non-trivial modifications in far, non observable regions of the plateau. On 
the other hand, if the potential has a flat inflection point, r/ c i can be large at intermediate 
wavelengths, too small to lie in the CMB observable window but still of astrophysical interest. 
This could have important consequences in possible non-linear effects at those small scales, 
such as the formation of primordial black holes (PBHs) [90]. In such models, the production 
of PBHs is calculated making use of the standard classical formulas for the amount of scalar 
perturbations. However, we have shown that in such regimes, stochastic effects largely modify 

19 Since v can only decrease, v < 10 _1 ° for all observable modes. 
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its value. An important question is therefore how this changes the production of PBHs in 
these models. In particular, it is interesting to notice that if the potential is concave (v" < 0), 
which is the case favored by observations [91-93], the leading correction in Eq. (4.4) is an 
enhancement of the power spectrum amplitude. However, as can be seen in the right panel 
of Fig. 4, as soon as one leaves the perturbative regime, this is replaced by the opposite 
trend: at the flat point, the classical result accounts for a diverging power spectrum while 
the stochastic effects make it finite. If generic, this effect may be important for the calculation 
of PBHs formation, and we plan to address this issue in a future publication. 

5 Conclusion 

Let us now summarize our main results. Making use of the SN formalism, we have shown how 
curvature perturbations can be related to fluctuations in the realized amount of inflationary 
e-folds in stochastic inflation trajectories. We have then applied “first passage time analysis” 
techniques to derive all the statistical moments of the number of e-folds, hence all scalar 
correlation functions in stochastic inflation. 

We have shown that the standard results can be recovered as saddle-point limits of the 
full expressions. The situation is therefore analogous to, e.g., path integral calculations. A 
new simple classicality criterion has been derived, which should replace the common estimate 
based on the ratio between the mean quantum kick and the classical drift during one e- 
fold. It shows that quantum corrections to inflationary observables are Planck suppressed in 
general (that is to say, they are proportional to V/M ^), but can be large if the potential is 
flat enough, even at sub-Planckian scales. For simple inflationary models where |dV/d<j!>|/U 
increases monotonously as inflation proceeds, the corresponding effects play a non-trivial 
role only at extremely large, non-observable scales. However, models containing a flat point 
in the potential between the Hubble exit location of the modes currently observed in the 
CMB and the end of inflation behave differently. First, the stochastic effects change the 
mean total number of inflationary e-folds and can therefore largely modify the location of 
the observational window along the inflationary potential. Second, the amount of scalar 
perturbations produced around the flat point is strongly modified by stochastic effects. This 
may be crucially important for a number of non-linear effects computed at these small scales, 
such as the formation of PBHs, or non-Gaussianity. 

Together with the case of tensor perturbations, which we have not addressed in this 
paper, we plan to study these issues in future publications. 
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A Why must we use the number of e-folds in the Langevin equation? 

In section 2, we have made explicit that different choices of time variable in the Langevin 
equation account for different stochastic processes. In this appendix, we explain why the cor¬ 
rect time variable to work with is the number of e-folds N , elaborating on already existing 
results. We first present a generic argument, based on a perturbation of the background equa¬ 
tions, before turning to an explicit comparison of stochastic inflation and QFT predictions, 
in order to identify the correct time variable. 


A.l Perturbations Equation derived from the Background Equation 


A heuristic derivation [17] of the Langevin equation relies on splitting the full quantum 
inflaton field into a coarse grained, large scales part tp, and a short-wavelength component 
0>, and on performing an expansion of the equation of motion in cj) > . As suggested in 
Ref. [37], the correct time variable should therefore be the one such that the equations for 
the perturbations, which must be established at the action level, can correctly be obtained 
from varying the equation of motion for the background itself, when written in terms of this 
time variable. In this section, we establish that this condition selects out N as the time 
variable. 

In the case where inflation is driven by a single scalar held </>, the action we start from 
is given by 


5 



'Ml 

2 


R - - v (</>) 


(A.l) 


From this action (and this action only), we first want to derive equations of motion for the 
scalar perturbations that can be compared with what will be obtained below from varying 
the background equation of motion itself. To make our point even more convincing, we go up 
to second order in the perturbations. This is why we expand the background fields {4>, g^p] 
at second order in the scalar perturbations. 20 When the time variable in the metric is the 
conformal time rj, one has 


<t> ( r\ , x ) = q i>(°) {rj) + q ( 77 , x ) + (??, x ) 


5oo — o' 


— 1 — 2a^ — 


9i0 = ~a 


diB {1) + -diB (2) 


(A.2) 


9ij = a 2 <j 5ij 1 - 2^ (1) - V> (2) 


+ 2 did 




1 


eW + -E (2) 
2 


The degrees of freedom introduced above are partially redundant and in absence of anisotropic 
stress, the scalar sector can be described in terms of a single gauge invariant variable. One 
possible choice is the Mukhanov-Sasaki variable [7, 42, 94] v, which can be defined, order by 
order, as the scalar held fluctuation 1 j)^ on uniform curvature hypersurfaces [95]. To hrst 

20 In this discussion, vector and tensor perturbations are irrelevant, which is why they are not taken into 
account. 
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and second orders, after a lengthy but straightforward calculation, one obtains [96] 
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(A.3) 
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2 [(#»>)"+ 




(A.4) 


From varying the expanded action, one can derive an equation of motion for the 
scalar perturbations, and in particular for a gauge invariant combination of them, say the 
Mukhanov-Sasaki variable. In this section, we want to compare this action based equation 
of motion for v with an equation of motion for the perturbation in (p coming from vary¬ 
ing the background Klein-Gordon equation. It is therefore important to work in a gauge 
where these two quantities, v and the perturbation in cp, are identical. By definition of the 
Mukhanov-Sasaki variable, this is the case in the uniform curvature gauge, where 


v (n) = 0 (n) 


(A.5) 


to all orders. In this gauge, one notably has ip = 0 to all orders [from Eqs. (A.3), (A.4) 
and (A.5), it is clear that this is at least the case up to second order]. 

The equation of motion for the scalar perturbations <p^ and <p ^ is therefore given by 
the one for t/ 1 ) and in this gauge. At leading order in the slow-roll approximation, and 
in the long-wavelength limit, they read 21 




V ' V " + \ (1) 2 
H 2 M 2 , 9H 4 M 4 t J 9 


(A- 6 ) 
(A.7) 


We now need to compare these equations with the ones that arise when varying the 
equation of motion for the background, and to find for which time variable they match. 


If t is used 

When cosmic time t is used, the leading order of the slow-roll approximation for the Klein- 
Gordon equation for the background is given by 


dcp _ V' 
d t ~ ”3 H(4>) ’ 


(A- 8 ) 


where we take H 2 ~ V/ (3M|j) at leading order in slow roll. When plugging <p = (f>W + cp^ + 
in this equation, one obtains at first and second order in the perturbations 




TTL + y ' 3 ) ^ 

2H 2 M 2 l l‘2H 4 Mp } J v 


(A.9) 

(A.10) 


21 In spite of the complexity of the field equations at second order, see e.g. Ref. [97], in the long-wavelength 
limit, it is sufficient [59] to use the local conservation of energy-momentum to establish Eqs. (A.6) and (A.7). 
Because this is not the main subject of this discussion, the corresponding calculations are not reproduced here 
but they can be found in Refs. [59, 95]. 
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One should stress that these equations do not apply to (/)(+ and since they are different 
from Eqs. (A.6) and (A. 7), which is why we use the notation />( 1,2 ) instead of (f>( 1,2 \ The 
differences with Eqs. (A.6) and (A.7) are displayed in red. One can see that several factors do 
not match. This is because in general, the equations for the perturbations must be derived 
from the action itself and cannot be obtained by simply varying the equation of motion for 
the background. 

If ds = H p a q dt is used 

For this reason, let us look for a time variable s which is such that the equations for the 
perturbations arise from varying the equation of motion of the background when written in 
terms of s. Let us assume that s is related to t thanks to a relation of the form 


ds = H p ( 4 >) a q (</>) dt, 


(A-11) 


where p and q are power indices that we try to determine. For example, when p = 0 and 
q = 0, s is the cosmic time t, when p = 1 and q = 0, s is the number of e-folds N, while 
when p = 0 and q = —1, s is the conformal time p. In terms of s, the equation of motion for 
the background is given by 


d(j> _ V' 
ds ~ ~3 HP+ 1 (</>) ’ 


(A.12) 


where again we take H 2 ~ V/(3M 2 l ) at leading order in slow roll. When plugging in 0 = 
0(°) + , one obtains at first and second order in the perturbations 
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H a 

1 ^ 



(A.14) 


Again, these equations do not apply to 0*- 1 - 1 and <p^ 2 ' > in general, since the correct ones are given 
by Eqs. (A.6) and (A.7) which is why we use the notation 0( 1,2 ). The differences between 
these two sets of equations are displayed in red. In order for the above to match Eqs. (A.6) 
and (A. 7), one must have q = 0 and (p+ l)/6 = 1/3, which gives p = 1, (p + l)/2 = 1, which 
also gives p = 1, and (p+ l)(p + 3)/36 = 2/9, which gives p = 1 or p = —5. As a conclusion, 
with p = 1 and q = 0 only, the equations for the perturbations (from what is shown here, up 
to second order in perturbation theory) can be seen as if they were derived from varying the 
equation of motion for the background. This choice corresponds to the number of e-folds N. 


A.2 Stochastic Inflation and Quantum Field Theory on Curved Space-Times 

To go beyond this generic argument, one can explicitly show [26, 27] that N is the time 
variable which allows one to consistently connect stochastic inflation with results from QFT 
on curved space-times. For example, let us consider the leading order of the fluctuations 
5(f> = p — cf >ci in the coarse grained inflaton field about its classical background value <j> c \. By 
“classical”, recall that we mean that 4> c \ is the solution of the equation of motion without the 


- 24 - 











noise term. We want to compute the mean square value of 5(f) and compare what we obtain 
with results coming from QFT calculations. For example, in Ref. [98], with a renormalization 
obtained by employing the adiabatic subtraction prescription on inflationary backgrounds, 
it is shown that in quadratic inflation where V = m 2 (f> 2 / 2, if 5(f) = 0 at time to, one has at 
leading order {see Eq. (48) of Ref. [98]} 

/ \ R 6 — ff 6 

= (A.15) 

where H means H ((f> c \) and Hq means H evaluated at time to- In the same manner, in 
Ref. [99], it was shown that in power-law inflation where a(t ) oc t p with p > 1, the same 
quantity is given by {see Eq. (29) of Ref. [99]} 

(W> - M 2 ) = ^ ( H o ~ H 2 ) . (A.16) 

Let us see how these results can be derived in the stochastic inflationary framework. We 
start from the Langevin equation (1.1) that we write 


d ip 

d N 




(A.17) 


where we have used H 2 ~ R/(3A l 2 ^) and where we recall that a prime denotes a derivative 
with respect to the inflaton field. Since (f> c \ is the solution of the above equation without 
the noise term, the noise term can be considered as a perturbation captured in 5(f). After 
expanding Eq. (A.17) in powers of 5(f), one gets for the leading order 5(f>^ 


d 5<f>W 
d N 





(A.18) 


Multiplying this equation by 5(f>^ and taking the stochastic average (•) leads to 


dN 


+ 4M P 2 , 




(A.19) 


In order to obtain a differential equation for (5(|)^ 1 ' ,2 } only, one needs to evaluate the right 
hand side of the previous equation. This can be done as follows. Letting = 0 at time 
Ao, a formal solution of Eq. (A.18) is given by 
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From this expression, since (£ (IV) £ (N 1 )) = 5 (N — N'), it is straightforward to see that 22 
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(w (,) ) = s ■ 


(A.21) 


This is why one obtains 


d N 


+ 4Mp2 ‘(§) (^ 1)2 ) 


H 2 
47r 2 


(A.22) 


2 “The 1/2 factor comes from the relation ffff f(x)S(x — X 2 )dx = f(x 2 )/%, which applies when the Dirac 

function is centered at a boundary of the integral. 
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Since the equation of motion for 4> c \ is simply given by dN = — H/ (2H'M 2 ^)d(j) c i, one can 
change the time variable from N to (j) c \ and the formal solution of the above equation can be 
written as 

(A ' 23) 

For quadratic inflation where H = m<p / (%/6M P1 ), this exactly gives rise to Eq. (A.15) while 
for power-law inflation where 23 H = Hq exp [— l/y/2p((/) — ^o)/Af P1 ], one exactly obtains 
Eq. (A. 16). Therefore, stochastic and standard field-theoretical approaches to inflation pro¬ 
duce the same results for the amount of field fluctuations. 24 

To emphasize the specificity of N as a preferred time variable, let us repeat the same 
procedure using the Langevin equation written in terms of t, 

(\th rx3/2 

jf = -2 MlH’ + —at) ■ (A.24) 

Since this corresponds to a different stochastic process as the one written in terms of N, we 
use again the notation 0 instead of (j>. At leading order in the noise, one obtains for 

d ^ + 2 M*H"8$n = (t) . (A.25) 

Again, multiplying this equation by 5(f>^ and taking the stochastic average leads to 

d/<h/>W 2 \ , 9 , tt 3/2 , , 

+ 4A/ 2 H” } = —— (t) <5<^ (1) ) . (A.26) 

In the same manner as before, making use of the formal solution to Eq. (A.25), 
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one can show that (t) ) = H 3 ^ 2 so that one needs to solve 



+ 4M 2 H" (^« 2 ) 


H 3 
47r 2 


Making use of the classical trajectory d t = —di?i c i/(2M| 1 i7 / ), one obtains 25 


(A.28) 



(A.29) 


23 One can show [100] that the potential associated with power-law inflation, for which a{t ) oc t p , is given 
by V(<j>) oc Since H 2 = V{</>)/(3M^) at leading order in slow roll, one obtains the given H(rj>) 

profile. 

24 Here we have established this result at leading order in perturbation theory. However, as shown in 
Refs. [17, 25], the stochastic approach can reproduce QFT results for any finite number of scalar loops and 
even beyond. 

25 This equation (A.29) also matches Eq. (13) of Ref. [101] where perturbative solutions of stochastic inflation 
are derived when formulated in terms of the cosmic time. 
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which is clearly different from Eq. (A.23). 2<) For example, for quadratic inflation, it reduces to 
) = 3(i^o — T7 4 )/(167r 2 m 2 ) which does not coincide with Eq. (A.15) and for power-law 
inflation, it reads (Scf)^ 2 ) = pH 2 / (Air 2 ) Ih(Hq/H), which does not coincide with Eq. (A.16). 

Finally and in passing, let us derive the corresponding results for the leading order of 
the mean fluctuation {5<j>). Since from Eq. (A.18) it is clear that ( dift 1 ') = 0, one has to work 
out ( 5q Expanding (p = c/) c \ + Seft 1 ^ + in Eq. (A.17), one obtains 

^ + 2 Ml (^pj ' 6^ + Ml (^pj " 6^ 2 = (N) • (A.30) 


When taking the stochastic average of the above equation, (5(j)( 1 ' ,2 } is given by Eq. (A.23) 
and (S(/>^£,) is given by Eq. (A.21), so that one obtains 


d(d^) 

dN 



(A.31) 


Using the classical trajectory d(t> c \ = —2MlH'/HdN, this equation can be written in terms 
of cj> c i, and after integration by parts, this gives rise to 
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where (d^ 1 ^) is given by Eq. (A.21). For example, when applied to quadratic inflation where 
V = m 2 (j) 2 /2, one obtains 
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which matches Eq. (49) of Ref. [26]. However, it is again worth noting that one would have 
obtained a completely different result starting from the Langevin equation written in terms 
of cosmic time t, namely 27 
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(A.34) 


This obviously differs from Eq. (A.32). 

To summarize the discussion, different time variables in the Langevin equation lead to 
different stochastic processes, and the only time variable which allows the stochastic inflation 
formalism to reproduce QFT calculations is the number of e-folds N. One should therefore 
always work with N when dealing with stochastic inflation. 
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